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At present, two different design concepts are being considered for 
the purpose of measuring heat flow as part of a Mars penetrator mission. 
The first of the tentative designs utilises temperature sensors emplaced 
along the trailing umbilicus at regularly spaced intervals, no greater than 
lm, which v/ill bo thermally coupled to the adjacent regolith radiatively 
and possibly couveciively (open hole configuration) or conductively (collap- 
sed hole configuration). Depending upon the depth of penetration, some or 
all of the umbilicus sensors will he subject to thermal perturbations due to 
diurnal and seasonal temperature variations, dissipation of initial heat and 
heat deposited during emplacement, and transients from energy released 
by the 10 or 20 watt KTG contained within the penetrator body. 

The second of ihe heat flow designs considered requires the radial 
deployment of two or more low thermal mass temperature sensors outward 
from the penetrator body over a vertical (depth) range on the order of Lm. 
Transient effects will, include the dissipation of initial heat of the sensor 
and accompanying radial deployment, arm and the eventually dominant ther- 


mal pulse originating from the penetrator initial heat and RTG power supply. 
Seasonal soil temperature variations may also be important if penetration 


depths do not exceed 15 m. 
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It i,M the purpose of thin report to describe resultH of numerical 
tiimul.'iUuiui used to iwveMtigate the feasibility and relative merits of each 
of the proposed design concepts tor a Mars penetrator heat flow measure- 
ment. It is to be noted that all discussion is of a theoretical and interpre- 
tutional nature. A complimentary investigation of the engineering and 
design problems associated with each of the proposed heat flow measure- 
ment modes will be essential and perhaps overriding in any final choice 
of experiment design. 


I. Umbilicus Sen sort! 

The effects of annual variations and RTG-induced thermal transients 
upon umbilicus temperatures in the open hole configuration have been pre- 
viously investigated (FY 1976 Progress Report on a Feasibility Study Eva- 
luating the Use of Surface Fenetrators for Planetary Exploration, NASA 
TM X-73, 181), Results of numerical models incorporating radiative and 
conductive thermal linkages indicated that the large annual subsurface 
temperature variations expected on Mars should be well coupled with the 
umbilicus sensors and could be accounted for given a sufficient time window 
of data (at least 0. 5 martian years), Significant RTG transient effects, 
propagated radiatively, will dominate umbilicus temperature histories 
within 1-2 m of the penetrator tail before one half martian year has elapsed. 
Thus, if annual variations are important at these depths, reliable mean 
temperature estimates will not be feasible. At depths where annual varia- 
tions are negligible (7 5-8 m), the necessary time frame of umbilicus 
data free of renetra<or transient effects would be reduced to the equilibration 
time of initial hea'-. deposition (on the order of tens of days), hi such a case, 
mearr'^giul mean temperature estimates could feasibly be made as near as 
lm from the penetrator tail. 
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Thts object of the present study of umbilicus models was to assesB 
the near surface effects of diurnal and annual temperature variations. Two 
specific effects were investigated: (1) The effect of a temperature dependent 
regolith thermal conductivity on the mean temperature profile; and, (2) The 
effect of diurnal variations propagating radiatively down an open hole. 

-- Temperature Dependent Thermal Conductivity - 

At the Apollo 19 and 17 lunar heat-flow sites, 35-45°K mean tem- 
perature differences were observed between the surface and depths greater 
than 10 cm. The interior heat flux can contribute only a negligible propor- 
tion of such gradients. Tile phenomenon can only be explained by the recti- 
fying effect of a strongly temperature dependent thermal conductivity in the 
upper few centimeters. During the warm lunar day, heat is transferred 
more effectively into the lunar surface than it can be removed during the 
cold lunar night. To conserve net flux over a lunation, a mean temperature 
gradient is established, confined mainly to the upper few centimeters. At 
the Apollo sites, variations of near surface thermal conductivity by factors 
of Z. 9-3. 3 were required to produce the observed mean temperature gra- 
dient. For the Mars situation, thermal conductivity temperature variation 
is expected to be limited by the dominance of the interstitial gas in the heat 
transfer process. Laboratory measurements, under simulated lunar and 
martian conditions (Fountain, J. A. and West, E. A. , 1970, J, G. R . 75, 
'>063-4069), resulted in a doubling of thermal conductivity under vacuum 
condition between 190 and 300 °K, but only a 30% increase over the same 
range in the presence of CC% at 7 mb. As an upper limit to the temperature 
dependence effect, we have assumed that the Mars regolith thermal conduc- 
tivity increases -v 80% between 160 and 290°K, a feasible surface temperature 
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varintion at near equatorial latitudes'. This corresponds to a variation of 
-V 2b% between 170 and 230 Q K, a feasible range of seasonal temperature 
variation in the upper meter. 

Mean temperature effects due to the diurnal heating cycle were 
calculated for two different levels! of conductivity. Results are shown in 
Figure 1. As; expected, over the short Mars diurnal period, (relative to 
the moon), the magnitude of the mean temperature change is small (3, 5 - 
6. fj°K) and limited to depths less than 10 cm. Even if a conductivity level 
ass high as 0. 001 w'cm°K characterized the upper layers, it is clear that 
martian diurnal mean temperature effects need not be considered as a per- 
turbing factor on steady-state gradient measurements; made below 20 cm. 

The much longer period annual variations on Mars, although of 
smaller amplitude, will penetrate to meter depths, and could result in 
significant mean temperature gradients as deep as two meters. The 
results shown in Figure 2 assume a seasonal surface temperature variation 
of - 200 1- 30 sin (n't) where « is the orbital frequency. Note that although 
the total effect is less than 1 °K, it extends over a depth range of v 2m, 
and is comparable to feasible steady- state gradients which could be produced 
by internal heat loss. Although the degree of conductivity temperature 
dependence employed in the models is most likely too high, the models do 
show that mean temperature gradients measured at depths shallower than 
two meters will require cautious interpretation. If umbilicus sensors 
were closely spaced ( v 0. 5 m), the degree of linearity of the mean tem- 
perature profile could provide one possible check on the reliability of 
interpreting a shallow mean temperature gradient in terms of internal heat 
lo as. 


- Open Hole Radiative Propagation of Diurnal Variations - 
Detailed modeling of the lunar heat flow probe-borestem geometry 



, '.bowed til.'? t radiative transfer of high frequency divrnal variations resulted 
hi enhanced propagation (relative to the adjacent regolith) of the diurnal wave 
along the open hole. The degree of amplification was proportional to the 
insulating proportion of the surrounding regolith, j, e, larger and deeper 
diurnal variations were produced for lower conductivity surroundings, In 
addition, small mean temperature differences were maintained between the 
Oorestem and adjacent regol’th to depths of 6 rt cm due to the lower lunation 
mean temperature of the hole "cap", i. e, the section of borestem which 
protruded above the lunar surface. 

Similar effects can be anticipated for a Mars penetrator open hole 

deployment with the afterbody analogous to the protruding lunar borestem. 

figure 3 illustrates the effect for a “worst case" model with a 9 cm diameter 

hole surrounded by a regolith of the lowest expected Mars thermal conduc- 
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tivity of k -- 2 x 10 \v/cm°K. The afterbody was modeled simply as a solid 
cylindrical section of perfect conductor capping the open hole with a height 
Id cm above the surface and a mass of 2, 3 kg. Results indicate that per- 
ceptible diurnal effects could be expected along the hole wall to depths of 
50 cm, even though diurnal variations in the undisturbed regolith are 
negligible at a depth of 20 cm. Such an effect must be anticipated if the 
hole remains open, and even time averaged borehole temperatures may 
not accurately reflect the regolith gradient at depths shallower than 50 cm. 

Radially Deployed Sensors 

Assuming sufficient radial extension ( 9 or more cm from penetrator 
body) and the limitation of impact-generated shock and deformation effects 
to within 1-2 cm of the penetrator, the primary factor affecting the measu- 
rement of preimpact regolith temperatures will be the thermal linkage from 
the penetrator body along the radial deployment arm (hereafter referred to 
as the “probe”) to the sensor. The penetrator and radial probe can be 


tlEPROJ)UCIBILTrY OF THE 

3©JK am * page is po- ■' 



- 6 - 


oxpeeted to maintain proimpact temperatures on the order of 50-1 50 °K 
above martian regolith temperatures, Impact generated heat will add to 
this difference, although most of the generated heat can be expected to 
reside in the penetrutor none and along a thin soil contact layer within 1mm 
of the penetrator body. The presence of a 10 or 20 W power source will 
prevent the penetrate r from equilibrating substantially to the regolith 
temperature. The primary objective of the models to be discussed is the 
determination of a feasible radial extension design which will allow an 
accurate measurement of the undisturbed regolith temperature gradient to 
be made before the penetrator thermal pulse significantly affects the 
preimpact thermal regime. 

Precise modeling of a cylindrical penetrutor with one or more radial 

extension arms requires a three dimensional geometry. However, since 

the penetrator -induced disturbance to the radial probe equilibration will be 

transmitted primarily u’ong the extension arm itself, the problem can be 

converted to two dimensions without significantly altering the early time 

history of penetrator-radial probe heat exchange. The models presented 

assume that the penetrator can be represented as an infinite heat reservoir 

at an elevated temperature V one hundred degrees above the undisturbed 

o 

regolith. The "infinite reservoir” assumption is justified by the much 
larger thermal mass of the penetrator body compared to the radial probe. 
The 100 °K temperature difference is assumed to include all preimpact 
and impact generated heat appropriate to the base of the radial probe. 
Schematics of the actual geometry and equivalent modeling geometry are 




Initial probe temperature is also as mimed to bo T , 100 °K above the 

surrounding regolith. The probe wa« assumed to be made of stainless steel 

^ 3 

with thermal condm tiv'fy k - 0. 13 ti w/c*m K and pc 1, 6wsec/cm °K. 

P r P 

The clmhc of ataiuler.H oteel in a ’’worst" cane in the sense that the large 
conductivity will enhance effects from the penetrator reservoir and the 
large volumetric heat capacity will slow equilibration to the surrounding 
oil temperature for a given probe diameter. However, the design requi- 
rement of a probe materiel of sufficient strength to allow penetration into 
highly compacted regolith may necessitate a choice of probe material 
of comparable thermal properties. 

The effects of four parameters were examined in the model runs, 
radial probe length I. , probe diameter R, regolith conductivity k, and 

probe conductivity, t . In all models a regolith heat capacity of^ =1.0 

3 1 

v/ sec /cm “K was assumed. Initial ’ probe temperature T 200 °K and 

regolith temperature T 100°K were used for convenience. All heat 

tx> 

transfer was assumed conductive. Results are shown in Figures 4-12 
where probe temperatures are shown at various distances from the pene- 
trator reservoir versus inverse time for three feasible values of regolith 


conductivity, K 


_4 

and 10 x 10 w/cm'K. Each figure is distinguished 


by the probe dimension.’ length and radius. The "s" data represent the 
transient solution for an infinite cylinder of perfect conductor equilibrating 
with a low conductivity infinite medium, and serve as a guide for deter- 
mining the onset and magnitude of the penetrator thermal pulse. At Iona 
times ( kt/Ryoc ' J ) the theoretical equilibration curves can be represented 
by: 

R yocp 

T(t) ' T CJU <r o ~ T J * Hcf" (R ■ probt ‘ radius> 


Thus, at sufficiently long times, if penetrator and end effects are 
negligible, a linear extrapolation of the inverse time plots to t-J-oo should 
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uruvicU: :ui accurate estimate of T even if T is not precisely lenown. 

The results of Figures 4-8 attempt lo show what combinations of probe 
length (i.i ), ratlin s (R) nud sons or location will allow accurate extrapo- 
lation to the regoiith equilibrium temperature T , before the effects of 
the thermal pulse from the penetrator "reservoir ' maslc the measurement. 
Futures ') and 10 show the effect of a nonnniform initial temperature dis- 
tribution along the radial extension arm. Figures 11 and 12 examine the 
.advantages lo be gained if the radial probe could be made of a much less 
conductive material than stainless steel, such as wound fiberglass. 

Specific results can best be discussed by addressing the Figures firectly: 

Figure 4: (L - 9.0; It - 0. 0f>) 

P 

Note that for a regoiith conductivity of 0.0002, even a sensor 
located at. the probe end (8. 7f3) begins to feel the effects of the penetrator 
pulse before an accurate extrapolation of the equilibration history could be 
made. The situation improves considerably for larger regoiith conduc- 
tivities, although penetrator effects begin to be felt at approximately the 
same time ( aj 1000 sec at the probe end). However, the probe end 
equilibration is sufficiently linear for the high lc - 0, 001 case that an extra- 
polation of the 200- 1000 sec curve would provide a T estimate accurate to 

oo 

0. FK. Note also, hov/.-ver, that end effects lead to more rapid cooling at 
the 8. 75 cm position than the theoretical case while at 7. 0 cm (2m from 
probe tip) the equilibration curve matches the theoretical cooldown almost 
exact! y before the onset of penetrator effects. This suggests that a sensor 
location slightly removed from the prubo tip may simplify the extrapolation 
to soil equilibrium even though the reequilibration will be most rapid at 
the probe tip. 

Figure 5: L - 9.0 cm; R 0. 15 cm 
P 

Note that, even for the high regolith conductivity case, the cooldown 
curves have been significantly disturbed by t - 800 sec. This, is not suf- 
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ficient time for the ft. 15 cm radius probe to substantially equilibrate, and 
extrapolation results, even for the high k cane, would be no better than + 1°K. 
If two radially deployed sensors separated by 1 m were to accurately measure 
a steady state gradient, a minimum absolute temperature accuracy of at least 
0. 1°K would be required. In actuality, it is feasible that the Mars regolith 
gradient could be of the order of 0, l°K/m, in which case extrapolation 
errors would be comparable to the quantity we are attempting to measure, 

i*'b>;u re (>: J, -- 20.0 cm ; R 0.05 cm 
P 

Here would be the ideal design for a radially deployed sensor. It is 

to be noted, of course, that a radial extension greater than the penetrator 

diameter (9cm) would severely complicate the engineering design in that 

the radial, arm would have; to be coiled or hinged in some manner prior to 

deployment. However, if such design complications could be overcome, 

it can be noon that any sensor positioning between 15 and 20 cm along the 

extension arm would provide equilibration histories essentially undisturbed 

by the penetrator pulse out to at least 2000 sec, a sufficient time span for 

0. l‘ K accurate T estimates to be made. Increased accuracy would be 
no 

possible with extended lime data as long as the sensor location remained 
unaffected by the penetrator pulse. Ideally, only the sensor precision need 
limit the accuracy of the T determination. 

Figure 7: (L 20.0; R - 0. 15) 

P 

Figure B: (L - 20. 0; R 0. 50) 

If the necessary rigidity could not be maintained in a 20 cm extension 

arm of 0. 05 cm radius, the effects of increasing the arm cross section 

will quickly lead to enhanced penetrator effect and the subsequent limitation 

of obtainable accuracy. Increasing the radius to 0. 15 cm may still 

allow a satisfactory T extrapolation from a sensor located at least 15 cm 

oo 1 



from the penetrator, H'W^vcr, considerably more than 2000 sec of equi- 
libration will he required, with the eventual limitation of penetrator pulse 
effects. This is seen moot clearly for the II 0. 5 ease (Figure 8) in which 
no extension arm location has equilibrated to within 10 “K of T . Even at 
the probe tip, (19.75) the first effects of the penetrator thermal pulse can 
tie seen a:, the cooldown curve crosscut the theoretical case at 2000 sec. 

It is clear that any feasible radial extension arm design would require a 
minimum possible probe radius since it can be seen from the long time 
theoretical expression (p 7) that the relative equilibration rate will scale 
inversely with the square of the- probe radius. Additionally, the conductive 
linkage to the penetrator along the extension arm will scale directly with 
the square of the probe radius 

Figure 9: L -- 20.0 ; R 0,05 

P 

Figure 10; L - 20.0; R - 0. 15 

p 

In Figures 9 and 10 we examine the reliability of the T extrapolation 
technique when the initial radial extension arm temperature varies with 
distance from the penetrator. Such a situation could he imagined if fric- 
ijonal effects are important during the radial extension process. The 
models shown represent an extreme case in which a 50 °K initial tempera- 
ture difference exists between probe tip and probe base. It can be seen, that 
with a sufficient time range of data, the initial differences will be smoothed 
and at extension arm locations greater than 12-15 cm a linear extrapolation 
of the 1/t cooling curves from 2000 sec out will provide comparable accurac 
of T estimates as wo re obtained for the initially isothermal probe cases 

!X> 

(Figures 6 and 7). It is to be noted, however, that conductivity estimates 
made from the slopes of the longtime IT curves require knowledge of the 
initial tempo i aUtre difference accurate to a comparable degree of the 
ie sired k precision. 


Linear temperature distribution 
in probe initially 
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Figure 11; L 20.0; U 

P 

0.05; K 

P 

K / 10.0 

ot; 

0. 0138w/cm°K 

Figure 12; L 20.(1, 3 

p 

5 0. 15; K 

P 

K / 1 0 . 0 

ss 

0.0138 w/cm°K 


lt ha a been suggested (Final Report and Recommendations of the Ad Hoc 
Surface Penetrator Sc ience Gurm.iitlee, 1976) that a much less conductive 
material, such as wound fiberglass, could be utilized for the radial probes. 
To examine the extension required for such a probe, equilibration curves 
have been plotted in Figures 11 and 12 fur various radial distances from 
the penetrator assuming a probe conductivity ten times less than that of 

stainless steel (K ), In Figure 11, a probe radius of 0.05 cm is assumed, 

s (■; 

and it can be seen that sufficient, undisturbed equilibration can be achieved 
ever the full range of expected Mars regolith conductivities by sensors 
located within 7-9 cm of the penetrator. An extrapolation accuracy of jh 
0, 01°K could ho achieved after one hour. This is to be contrasted with the 
stainless steel probe case (Figure 6) in which sensor placements = 15 cm 
would be required for the lower range of expected regolith conductivities. 

The advantages of requiring a sensor extension no greater than the pene- 
trator diameter are obvious from the design viewpoint. 

Figure 12 illustrates the low probe conductivity case for a probe 
radius of 0. 15 cm. Again, the delay of penetrator influence can be seen 
when compared to the comparable stainless steel case of Figure 7. Sensors 
located at the critical 9 cm position would equilibrate free of penetrator 
effects for periods greater than one hour, allowing extrapolation accuracy 
to an i" fimated i 0. l r 'K, which would be adequate for all but the lowest 
expected values of Mars steady state regolith gradients. 

Tri summary, the analysis of purely conductive effects indicates that 
the possibility of measuring the Mars regolith temperature gradient by the 
radial extension method is constrained primarily by design limitations. 
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To jU'M.'.un? undisturbed regolith temperatures to _+ 0. 1°K will require 
the deployment of a rarlial extension arm with a high length/coos section 
ratio . If sensor extension is limited to 9 cm (penotrator diameter), a low 
conductivity probe, such as wound fiberglass, would yield extrapolation 
iiecuracy of 0. 01 °K for a 0, 05 cm arm radius and + 0. 1 K for a 0, 15 cm arm 
radius. If a more conductive material were used, radial arm extension 
greater than 9 cm would be required to insure comparable extrapolation 
accuracy over the full range of expected regolith conductivities. If the 
necessary design could be implemented, the success of a heat-flow mea- 
surement would then depend on the sufficiency of penetrator emplacement 
depth. If significant annual variations are present at the measurement 
depth, it will not be possible to distinguish periodic fluxes from the flux 
from the interior because of the very small time allowance for unperturbed 
measurement. The situation worsens for the higher feasible regolith con- 
ductivity values, an shown in Figure 13. For a given heat flow, the steady 
state gradient will decrease inversely with the conductivity. For example, 

if Mars heat flow is comparable to the measured lunar heat flow (2, 0 
2 

gw/cm ) transient gradients comparable to the steady state gradient will 
persist to depths of 7-9 m if the regolith conductivity is as high as 0.001 
w/oni^K. Such conductivity values have been inferred for selected martian 
regions from Mariner IR results. Even for the most favorable situation, 
k - 0.0002 w/cm°K, significant annual effects could persist to 4 m. It 
should be emphasized that even a theoretical correction for annual effects 
would not be possible without a reliable estimate of the thermal diffusivity 
profile from the surface to penetrator depth. 

Conclusions and Recommendations 
The relative merits of the umbilicus and radial extension techniques 
for the measurement of Mars heat flow may be summarized as follows: 
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Umbilicniti Deployment : 

1) .Sufficient time for analysis of transient and .seasonal elfects for 
f I'UiJorn at leant. 1 m removed from the penetrator tail. Effects duo to 
diurnal variations withit) tin* borehole will likely obneure interpretation of 
any measurements within 50cm of the surface. If regolith conductivity is 
sufficiently temperature dependent ( v 25°' variation over a 60°K range), 
a mean temperature gradient unrelated to interior heat flow could produce 
significant interpretation ambiguity to depths of 2 m, Sxich a dependence is 
not considered likely in view of the interstitial gac presence on Mars, 

2) I< 'or sufficient penetration ( .*■ 5 rr.) and sensor spacing -■ lm, 
three or more* mean temperature mes surnjn mV* could provide added 
confidence in the estimate of regolith gradient. Possible depth variations 
of thermal properties could be? detected. pruvirRrg an important constraint 
on any non-linearity j, sweat in the gradient measurement. 

3) Even in the absence of sufficient penetration to allow a reliable 
steady state gradient measurement, one or more umbilicus sensors could 
still provide data constraining regolith thermal properties and the diurnal 
and seasonal thermal regime. 

Radial Extension Deployment: 

1) Excellent thermal contact with the regolith material. If sufficient 
extension of a low conductivity radial probe could be achieved, a highly 
accurate measurement of regolith temperature, uncomplicated by radiative 
and convective thermal linkages, could be made. II should be pointed out 
that numerical models end convective simulation models (T. Canning, 

NASA Arnes) appear to demonstrate that umbilicus sensors, adequately 
removed from penetrator and surface effects, should be well coupled with 
the adjacent regolith. However, the conductive coupling of a slender 
probe within a particulate soil has been repeatedly demonstrated by labora- 
tory measurements, and must be considered superior in this respect. 
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Z) Inherently greater measurement depth for a given penetration. 

A situation can be imagined in which insufficient penetration precludes 
a reliable gradient measurement at umbilicus depths due to annual 
variations which are quickly masked by penetrator induced transients, 
but for which annual variations at the penetrator depth are negligible. In 
tiuch a case, only the radial extension design could provide a reliable 
: ; rad ient mea our ement . 

In summary, it is our feeling that the narrow time constraint of the 
radial extension design is the overriding negative factor in the choice of 
alternate heat flow measurement techniques. The lunar experience clearly 
demonstrated the value of a multisenaor network which could measure 
transient and periodic temperature variations as well as the steady state 
gradient. Unless new information is obtained demonstrating the likelihood 
of large penetration depths ( > 5 m) over much of the martian surface, no 
design concept which has no chance of measuring seasonal variations can 
be considered favorable. 

The ideal experiment, of course, would include both umbilicus and 
radially deployed temperature sensors. We would therefore still argue 
foi a serious engineering study of the design feasibilities and limitations 
of both concepts. 
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Figure Captions 

Figure 1 Diurnal mean temperature versus depth for two models o£ 

Mars regolith thermal conductivity temperature dependence. 
Heat flow, Q ~ 0, is assumed. 

Figure 2 Annual mean temperature versus depth for the two models of 

Mars conductivity used in Figure 1. The dashed line illustrates 
a feasible gradient that could be produced by the interior heat 
loss, 

Figure 3 Diurnal variations expected in a 9 cm diameter open hole on 

Mars compared with fluctuations at comparable regolith depths. 

Figure 4-8 Inverse time equilibration curves for sensors deployed radially 

outward from the penetrator along a steel extension arm. Each 

figure is distinguished by the radial extension arm length (D ) 

P 

and radius (R). Numbers on each curve identify the radial dis- 
tance (in cm) along the extension arm from the penetrator body. 
In all cases, a preimpact regolith temperature of 100 °K was 
assumed with the r? dial probe initially at 200 °K and the pene- 
trator body maintained at 200 “K, Each model was run for the 
three feasible values of Mars regolith conductivity shown in the 
figures. 

Figures 9,10 Equilibration curves for the same models as shown in Figures 
6 and 7 with the exception that a uniform radial gradient of 
2, 5°K/cm exists initially along the radial extension arm. 
Numbers in parentheses identify the initial temperature at the 
hypothetical sensor location. 
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Figurea 1 


Figure 13 


1, 12 Equilibration curves for the same models as shown in 
Figures 6 and 7 with the exception that a probe conduc- 
tivity ten times less than stainless steel is assumed. 

Mars regolith temperature profiles at two different 
annual phases for two homogenous models of feasible re- 
golith diffusivities, A simple harmonic variation has been 
astmmed at the surface of 30°K amplitude. 
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